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1. INTRODUCTION 


This report is the final report for NASA Grant NAGW-821 entitled “Kinetics of Thermo- 
chemical Gas-Solid Reactions Important in the Venus Sulfur Cycle”. This report covers the grant 
period of 11/01/85 - 10/31/88. 

The major accomplishments achieved during this grant period are described in section (2a) of 
this report and in the enclosed reprints and preprints. 


2. FINAL REPORT 

(a) Accomplishments 

The thermochemical net reaction 

CaCO$ (calcite) + S0 2 (p) —* CaSO± (anhydrite) -\-CO(g ) (l) 

is predicted to be an important sink for incorporation of S0 2 into the Venus crust (Prinn 1985a, 
b, c; Fegley 1988; Fegley and Prinn 1988). Furthermore, several authors have noted that the C0 2 
partial pressure at the Venus surface ( Pco , ~ 92 bars) is almost identical to the equilibrium C0 2 
pressure (Pco, ~ 102 bars) developed at the Venus surface temperature of 737 K from the Urey 
equilibrium 


CaC0 3 (calcite) -+- St0 2 (quartz) = CaSiO$ (wollastonite) + C0 2 (<?) (2) 

and have proposed that this equilibrium is buffering the C0 2 pressure in the Venus atmosphere 
(e.g., Urey 1959; Mueller 1964; Lewis 1970; Vinogradov and Volkov 1971). These predictions 
concerning possible roles for calcite on Venus, the exemplary nature of reaction (1), and the 
availability of high purity CaC0 3 in severed forms led us to study reaction (1) as a paradigm 
for SO 2 removal on Venus. 

Our philosophy in studying this reaction was to establish the reaction rate law and to under- 
stand the dependence of rate on experimental variables such as temperature, SO 2 partial pressure, 
C0 2 partial pressure, and 0 2 partial pressure. We can then extrapolate the experimental results 
to conditions of pressure, temperature, and composition relevant to the surface of Venus. However, 
our objective was not to attempt a precise simulation of the Venus surface conditions. For exam- 
ple, although the SO 2 molecular number density at the surface of Venus was reproduced in our 
experiments, the C0 2 molecular number density (i.e., total pressure) was not. Instead, an inert 
dilutant gas (Ar) was used to study the variation of the rate with the C0 2 partial pressure. 

Our experimental approach was a variant of the thermogravimetric method employed by ma- 
terials scientists to study the kinetics of many thermochemical gas-solid reactions (e.g., Beruto and 
Searcy 1974; Bratton and Brindley 1965; Gordon and Kingery 1967; Powell and Searcy 1978; Turk- 
dogan 1968; Turkdogan et al., 1965). The principal difference was that we heated different samples 
for varying time periods along isotherms instead of continuously monitoring the weight change 
of a single sample as a function of time. Our method was preferable in this case because of the 
slow reaction rate (days to weeks) and because it allows us to use powerful microscopic techniques 
such as scanning electron microscopy (SEM) and energy dispersive spectroscopy (EDS/SEM) to 
examine the surface chemical and morphological changes as the reaction proceeds. Results of these 
techniques will be discussed and illustrated later. 
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Clear calcite crystals (Iceland spar) were heated at constant temperature in S0 2 -bearing 
gas streams for varying time periods. The gas-solid reaction apparatus, which is schematically 
illustrated in Figure A1 of the Appendix, was used to do this. The purity of the calcite crystals 
was checked by chemical analysis using inductively coupled plasma emission spectroscopy (ICP). 
The analytical results, which are given in Table A1 of the Appendix, show that total impurities 
are ~ 350 ppm by mass with Mg (~ 200 ppm) being the major impurity. X-ray diffraction was 
also used to check the calcite crystals. These results (Table A2) show that all the observed X-ray 
lines match the literature values for calcite. 

Before heating, the calcite crystals were carefully weighed using a Mettler analytical balance 
with a 0.01 milligram sensitivity. Their geometric dimensions were measured to the nearest 0.00025 
cm (0.0001 inch) with a micrometer. These measurements are necessary to calculate the surface 
areas of the samples because we are interested in measuring rates per unit area. The micrometer 
measurements were checked by comparing the observed sample weights with the geometrically 
calculated weight. These values generally agreed within ± 5%. Also, replicate weight and sample 
size measurements were made to calculate 1 a uncertainties via standard methods of small number 
statistics (Shoemaker et al., 1974). 

The experiments were done at ambient atmospheric pressures with analyzed S0 2 - C0 2 gas 
mixtures nominally containing 1% (by volume) S0 2 . These mixtures were chosen so that the S0 2 
molecular number density in the laboratory experiments is about the same as at the surface of 
Venus where the atmospheric pressure is about 100 times larger (Seiff 1983) but the S0 2 volume 
mixing ratio is about 100 times smaller (von Zahn et al., 1983). As Table A3 shows, the mixtures 
also typically contained < 1 ppm 0 2 and < 30 ppm H^O as impurities. Standard flow meters were 
used to monitor and control gas flow during the experiments. A microprocessor controlled vertical 
tube furnace (Deltech model DT-31-RS-OS) was used in all experiments. Separate thermocouples 
were used to control the furnace (Pt^Rh^ vs. PtjoRhzo) and to monitor the sample temperature 
( Pt vs. PtgoRhio). The sample thermocouple was calibrated at the melting points of Au and Ag 
prior to use and was placed a few mm away from the samples in the constant temperature zone (~ 
40 mm long) of the furnace. The combined uncertainties in temperature control and measurement 
are less than ± 5°C. At the end of an experimental run, samples were quenched in the furnace 
by simultaneously shutting off the power and switching from the S0 2 - C0 2 gas mixture to pure 
C0 2 . (A C0 2 pressure of 1 atmosphere will prevent calcite dissociation at temperatures less 
than 1175 K, which were used in our experiments.) There are several good reasons for using this 
quenching procedure: (a) the slow rate of reaction leads to negligible reaction during the ~ 10 
minutes required to flush the furnace tube with pure C0 2 , (b) the initial cooling rate of 25-30° C 
per minute, depending on the run temperature, and the strong dependence of rate on temperature 
also leads to negligible reaction during the — 3-4 minutes required to cool the sample by 100° C. 
We note that for our highest temperature run used to determine the activation energy for reaction 
(1), any additional reaction during the first 4 minutes of cooling, when substantial amounts of the 
reactant gas are still in the furnace tube, is less than 0.04% of the total number of S0 2 molecules 
reacted during the entire run (see experimental conditions for R43 in Table A5). 

The rate of reaction (1) was determined by three independent methods: (a) reweighing the 
reacted samples and measuring the weight gain at the end of reaction, (b) using SEM to measure 
the porosity and thickness of the reacted surface layers from micrographs of external and fracture 
surfaces, and (c) chemical analyses by ion chromatography of the amounts of SO^~ in the reacted 
samples. 

SEM examination of many samples (in combination with X-ray line scans and EDS/SEM) 
showed that CaSO± layers covered all external surfaces of the reacted calcite crystals. The X- 
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ray line scan across the fracture surface of a reacted sample (Figure A2) and the fracture surface 
micrograph of another reacted sample (Figure A3) are representative of the results obtained. X-ray 
diffraction patterns of reacted calcite crystals and of more extensively reacted CaCO 3 synthetic 
powders showed only lines due to anhydrite and to unreacted calcite. This is illustrated in Table 
A4 which shows that all (non-calcite) lines in the product are due to anhydrite. Both the literature 
values for anhydrite and our separate results for a natural anhydrite sample from the MIT mineral 
collection agree with the values for the product. 

The experimental rate data and run conditions for the subset of runs used to determine the 
activation energy for reaction (1) are summarized in Table A5 and Figure A4. (Other experimental 
runs were either preliminary runs, runs used to derive the rate dependence on other experimental 
variables, or were runs on other calcium minerals.) These rate data are the first experimental 
measurements of the rate of any chemical weathering reaction postulated to occur on 
the surface of Venus. The use of these data in estimating the present-day rate of volcanism on 
Venus will now be described. 

A weighted linear least squares fit to all the rate data in Table A5 gives the equation for the 
rate R = io 19 - 64 l ±0 - 28 ) exp (— 15, 248(±2970)/T) molecules cm -2 sec -1 . As illustrated in Figure 
A4, the corresponding activation energy is 126.8 ± 5.4 kJ/mole. The rate of reaction (1) was 
extrapolated downward to Venus surface temperatures (~ 660 — 750 K) using this equation, the 
Pioneer Venus radar altimetry data (Masursky et al., 1980), and the Pioneer Venus atmospheric 
(P, T) profile (Seiff, 1983) to take the altitude dependence of the rate into account. 

The results of these calculations Eire displayed in Figure A5, which illustrates several important 
points. First, because the rate of reaction (1) exponentially decreases with decreasing temperature, 
the reaction of S0 2 with calcite on the Venus surface will be slower at higher elevations (and thus 
lower temperatures) than it will be at lower elevations (and thus higher temperatures). Second, 
assuming that calcite is uniformly distributed over the Venus surface, most of the global anhydrite 
formation via reaction (1) will take place at and around the median planetary radius of 6051.5 km. 
This is simply a consequence of the unimodal surface topography distribution on Venus. Third, 
summing up the reaction rates in all elevation increments (see Fig. A5) gives a derived global 
mean S0 2 reaction rate ~ 4.6 x 10 10 molecules cm -2 s -1 . Based on the median planetary radius 
of 6051.5 km (Masursky et al., 1980), this corresponds to a SO 2 reaction rate of ~ 2.1 x 10 29 SO 2 
molecules s" 1 over the entire Venus surface. 

Chemical analyses of the Venus surface by the Venera 13, 14 and Vega 2 spacecraft (Surkov 
et al., 1984, 1986) gave a CaO content (by mass) of 7.90% (weighted mean). Assuming that 
our experimentally determined rate is representative of the SO 2 depletion rate from the Venus 
atmosphere via reaction with calcium minerals on the Venus surface, the observed S0 2 column 
density of 2.2 x 10 23 molecules cm -2 (von Zahn et al., 1983) would be removed from the Venus 
atmosphere in ~ 1.9 x 10 6 years in the absence of a comparable sulfur source. This is a geologically 
short time period. However, unless the observed clouds, which are formed by the sunlight- 

powered conversion of S0 2 into H 2 SO 4 cloud particles (Esposito et al., 1983), are equally short- 
lived, their maintenance requires a comparable sulfur source. 

This source must be endogenic because the flux of sulfur provided by infalling cosmic material 
is too small to balance the sulfur flux into the Venus crust. The calculated S0 2 removal rate 
over the entire surface of Venus is equivalent to ~ 2.8 x 10 13 grams of sulfur per year. However, 
the measured terrestrial flux of infalling material, which is roughly applicable to Venus because of 
the similar sizes and masses of the two planets, is only ~ 7.8 x 10 10 grams per year (Kyte and 
Wasson 1986), and only contributes at most ~ 4.6 x 10 9 grams sulfur per year assuming the most 
sulfur-rich Cl chondrite composition (Mason 1971). 
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The most plausible endogenic sulfur source is volcanism, which has occurred on Venus in the 
past (Basilevsky and Head 1988), and which may also have led to the increased S0 2 levels above 
the Venus cloud tops that were observed by the Pioneer Venus orbiter (Esposito 1984; Esposito et 
ai, 1988). The rate of volcanism required to balance the atmospheric SO 2 depletion (via reactions 
with calcium minerals on the Venus surface) depends on the sulfur content of the erupted material 
(both gas and magma). If this material has an overall S/Si mass ratio of 0.03, the weighted mean 
of the Venera 13, 14 and Vega 2 surface analyses (Surkov et al ., 1984, 1986), the corresponding 
rate of volcanism is ~ 1 km 3 yr -1 . (This model as well as the other two models below assumes 
a rock density of 3 g cm -3 .) Figure A6 illustrates that this is about 20 times smaller than the 
terrestrial value. 

Two other plausible models for the S/Si ratio of the erupted material on Venus are also shown 
in Figure A6. The first model assumes that the bulk S/Si ratio of the erupted material is ~ 0.1, as 
in ordinary chondrites (Mason, 1971). The corresponding volcanism rate is ~ 0.4 km 3 yr -1 . The 
second model assumes that the S/Si ratio of the erupted material is ~ 0.004, as in the Earth’s 
crust (Ronov and Yaroshevsky 1976) and gives a volcanism rate of ~ 11 km 3 yr -1 . Thus, three 
plausible models for the bulk S/Si ratio of the erupted material give volcanism rates ranging from 
~ 0.4 to ~ 11 km 3 yr -1 , a range that is smaller than the terrestrial rate of ~ 20 km 3 yr -1 . 

All three rates will scale as e -1 , where e is the S degassing efficiency (0 < e < 1) of the 
erupted material. This efficiency is probably close to unity because sulfur gases will be expelled 
directly into the atmosphere and the hot, sulfur-bearing magmas will rapidly react with and lose 
S to the atmosphere, both during the eruption and later during cooling. An example of the facile 
nature of these S degassing reactions is provided by our preliminary experiments on pyrite ( FeS 2 ) 
weathering. Figure A7 shows pyrite crystals on a natural dolomite ( CaMg(C0 3 ) 2 ) crystal prior 
to any high temperature reaction. However, as Figures A8 and A9 show, heating at ~ 980° C for 
only 3 hours is sufficient to degas the sulfur from the pyrite crystals and convert them to iron 
oxides. Reactions at typical magmatic temperatures of ~ 1100-1200° C will be even faster and 
sulfur degassing will also be even more thermodynamically favorable. 

Figure A6 also illustrates the two widely discordant volcanism rates that have been estimated 
by geophysical methods. Solomon and Head (1982) and Turcotte (1988) have both estimated 
volcanism rates in the range of ~ 200-300 km 3 yr -1 by scaling terrestrial heat production to Venus 
and then assuming that volcanism is the sole heat transport mechanism. However, Grimm and 
Solomon (1987) have estimated a much lower limit of less than 2 km 3 yr -1 for the Venus volcanism 
rate by an analysis of the areal density of impact craters on the Venera 15, 16 radar images. 

The geochemically derived volcanism rate, which is based on our experimental study of the 
SO 2 + calcite reaction, is independent of these two discordant rates. In fact, the geochemically 
derived rate can in principal be used to test if either (or neither) of the geophysically derived rates 
is correct. Finally, our work also suggests that Venus is less volcanically active than the Earth. 


(b) Publication Record 

The following refereed papers supported by NASA Grant NAGW-821 during the 3-year 
period 11/01/85 - 10/31/88 have been published or are in press: 

(1) R.G. Prinn and B. Fegley, Jr. “The Atmospheres of Venus, Earth, and Mars: A Critical 
Comparison,” invited chapter in Ann. Rev. Earth Planet. Sci., 15, 171-212 (1987). 

(2) R.G. Prinn and B. Fegley, Jr. “Solar Nebula Chemistry: Origin of Planetary, Satellite, and 
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Cometary Volatiles,” invited chapter in Planetary and Satellite Atmospheres; Origin and 
Evolution (University of Arizona Press, Tucson), in press for 1988. 

(3) B. Fegley, Jr. “Cosmochemical Trends of Volatiles in the Solar System,” invited chapter 
in Workshop on the Origins of Solar Systems ed. J.A. Nuth and P. Sylvester, LPI Technical 
Report No. 88-04. 

The following manuscripts supported by NASA Grant NAGW-821 have been submitted 
for publication or are in preparation: 

(4) B. Fegley, Jr. and R.G. Prinn “The Rate of Volcanism on Venus,,” submitted to Nature. 

(5) B. Fegley, Jr. “The Rate of Chemical Weathering on the Surface of Venus: Experimental 
Study of Anhydrite Formation,” in preparation for the H. Wanke Festschrift volume of Z. 
Naturfor. 

Also, the following papers based on research supported by NASA Grant NAGW-821 were 
delivered at scientific meetings (speaker underlined) during the past year: 

(6) B. Fegley. Jr. “First Laboratory Measurements of Thermochemical Kinetics of S0 2 Reactions 
with Putative Venus Crustal Minerals,” 19th Annual Division of Planetary Sciences of AAS 
meeting, Pasadena, CA, 10-13 November 1987. Abstract: Bull. Am. Astron. Soc., 19, 871. 

(7) B. Fegley. Jr. “Thermochemical Kinetics of S0 2 Reactions with Possible Venus Crustal Min- 
erals: First Data for Calcite,” 19th Lunar and Planetary Science Conference, Houston, TX, 
14-18 March 1988, Abstract: Lunar Planet. Sci., XIX, 315-316. 

(8) B. Fegley. Jr. “First Laboratory Data on the Rates of Sulfur Gas-Mineral Interactions Relevant 
to Venus,” Brown University- Vernadsky Institute Microsymposium 7, Providence, RI, 21-23 
March 1988. 

(9) B. Feglev. Jr. “The Rate of Chemical Weathering on the Surface of Venus: Experimental Study 
of Anhydrite Formation,” scheduled to be presented at the Wanke Festschrift Conference, 
Max-Planck-Institut fur Chemie, Mainz, W. Germany, 5-6 September 1988. 

In addition, the PI was invited to contribute a report based on research supported by NASA 
Grant NAGW-821 to the 1988 Planetary Geosciences Annual Report. This is the first time 
this report, which is intended to publicize the research results in Planetary Geosciences, has been 
prepared. Approximately 40 Pi’s were invited to contribute. 


(c) Students Involved in Research 

The following MIT undergraduate students participated in research supported by NASA 
Grant NAGW-821. The student’s contributions to specific publications from section 2b are 
listed in parenthesis after the student’s name. 

(1) M. Camitta (assisted in building gas-solid reaction apparatus) 

(2) A. Petropoulos (paper 4) 

(3) S. Rowell (paper 4) 

(4) B. Vanderlaan (paper 2) 
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TABLE A1 


Impurities in Calcite Crystals (ICP Analyses)* 
Impurity Level (ppm by mass) 


Element 

SamDle 1 

SamDle 2 

Al 

6.48 

6.92 

Mg 

220 

212 

Mn 

24.3 

22.1 

Si 

12.4 

11.8 

Sr 

72.2 

69.8 

Zr* 

23.6 

26.2 

OTAL IMPURITIES 

359 

349 


f Inductively Coupled Plasma Emission Spectroscopy 

t The Zr impurity probably comes from the ZrO^ grinding media used to pulverize the crystals 
for chemical analysis. 



table 


A2 


fieactaru 

3.025 

2.485 

2.280 

2.089 

1-905 

1.872 


X-R a y 


a ctj 0o 


Caicite 
d' Valu eg 


Results 


^eratyre 

3-035 

2.495 


2-285 

2.095 


1-913 
1.87 5 



TABLE A3 


Analyses of SO 2 - CO 2 Gas Mixtures® 


RUN 

% SOo 

vvm O ' > 

v vm HoO 

R36 6 

1.06 

< 1 

21 

R43 6 

1.06 

< 1 

21 

R44 6 

1.06 

< 1 

21 

R45 c 

1.08 

< 1 

34 

RSO'* 

0.96 

< 1 

11 

R51 d 

0.96 

< 1 

11 


0 Only the runs used to determine the activation energy of reaction (1) are listed. 

6 The same tank was used in these runs and the impurity analyses were done when the tank 
was nearly empty. 

c The impurity analyses were done when the tank was nearly empty. 
d The same tank was used in these runs. 



TABLE A4 


X-Ray Diffraction Results 

Anhydrite 

d-Values 


Product 

Anhydrite 

Literature 

3.484 

3.479 

3.498 

2.835 

2.837 

2.849 

2.318 

2.320 

2.328 

2.198 

2.203 

2.208 

2.173 


2.183 

1.985 


1.993 

1.742 

1.745 

1.749 

1.643 

1.643 

1.648 


• All lines in product can be attributed to anhydrite or unreacted calcite 

• 23 hour exposure on product, 12 hour exposure on natural anhydrite 



TABLE A5 


Rate Data Used to Calculate the Activation Energy 


Sample 

T(°C) a 

Duration 

(Hours) 

Rate x 10 -12 
( molecules cm -2 : 

R44-1 

602 

95 

1.67 ± 0.19 c 

R45-1 

602 

314 

0.85 ± 0.77 

R45-3 

602 

314 

1.05 ± 0.12 c 
1.01 ± 0.11 c 

R45-4 

602 

314 

1.20 ± 0.51 

R36-1 

746 

96 

0.68 ± 0.48 d 
8.43 ± 2.93 

R36-2 

746 

96 

13.6 ± 6.0 d 
23.2 ± 6.2 

R50-1 

747 

192 

19.9 ± 11. 0 d 
9.19 ± 1.60 

R50-3 

747 

192 

10.5 ± 7.8 

R51-1 

748 

96 

13.2 ± 4.0 

R51-2 

748 

96 

10.3 ± 5.6 

R43-1 

848 

91 

60.4 ± 6.8 

R43-2 

848 

91 

57.1 ± 31. 7 d 

53.2 ± 7.0 

R43-3 

848 

91 

67.2 ± 43.4 d 
58.1 ± 7.0 

Notes to Table 1 



48.4 ± 12.1 d 


a The estimated temperature uncertainty is ± 5°C. 

6 The number of SO 2 molecules reacted per cm 2 surface area per second. The rate is from 
weight gain measurements unless noted otherwise. The 1 a uncertainties are also listed. 

c The rate from ion chromatography analyses of SO\~ content of reacted sample. 
d The rate from scanning electron microscope (SEM) measurements of anhydrite layer thickness 
and porosity. 
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Figure A2. Scanning electron microscope (SEM) micrograph of the fracture surface of a reacted 
calcite ( CaCO $ ) crystal. The scale bar is 50 micrometers long. The horizontal white line on the 
micrograph shows the position of an X-ray line scam for the element sulfur. The way white line 
shows that sulfur X-rays are produced only at the reacted surface where grains of the mineral 
anhydrite ( CaSO 4 ) are formed as a result of the gas-solid reaction. 
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Figure A3. A SEM micrograph of the fracture surface of another reacted calcite crystal. The 
scale bar is 100 micrometers long. All external surfaces of the crystal are covered with a layer of 
anhydrite ( CaSO 4 ) grains. This layer gradually becomes thicker as the gas-solid reaction continues. 
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Figure A4. Arrhenius plot of rate data from Table A5 vs. 1/T. The line is a weighted linear least 
squares fit. Rates determined from independent methods (wt. gain, SEM, SO+~ analyses) for the 
same sample agree within the 1<7 experimental uncertainties. 
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Figure A5. The rate of reaction (1) as a function of temperature and elevation on Venus. The 
bar graph shows the % of the total rate in each elevation increment. The solid line shows how the 
rate decreases with decreasing temperature. The global mean rate is ~ 4.6 x 10 10 molecules cm -2 
sec -1 
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Figure A6. Independent geochemical and geophysical estimates of the rate of volcanism on 
Venus are compared to the terrestrial rate. Three plausible models for the S/Si ratio of erupted 
material on Venus yield rates of ~ 0.4 to ~ 11 km 3 yr _1 . Adopting the Venus surface S/Si ratio 
from Venera 13, 14 and Vega 2 analyses yields a rate of ~ 1 km 3 yr -1 , or about 20 times less than 
the rate of volcanism on Earth. 
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Figure A7. Scanning electron micrograph of pyrite (FeSz) crystals on a dolomite CaMg(C0 3 ) 2 
substrate. The scale bar is 100 microns long. See Figure A8 which shows pyrite crystals after high 
temperature reaction with CO 2 (gas). 
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Figure A8. Scanning electron micrograph of reacted pyrite ( FeS 2 ) crystals heated at ~ 980° C for 
3 hours. The scale bar is 100 microns long. The “fluffy’ particles are the reacted pyrite crystals 
which have lost all sulfur and have been oxidized by the CO 2 atmosphere in the furnace. Compare 
to the unreacted crystals in Figure A7. 
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Figure A9. Chemical analysis of the reacted pyrite crystals in Figure A8. The analysis is done by 
energy dispersive spectroscopy on the scanning electron microscope (EDS/SEM). Peak intensities 
are on the vertical axis (1000 counts full scale) and peak energies are on the horizontal axis (0- 
10 Kev). Note that the reacted pyrite crystals have lost all sulfur and are now Fe oxide. Only 
Fe peaks from the reacted pyrite and Ca and Mg contamination from the surrounding dolomite 
(partially converted to CaO + MgO ) show up because the oxygen and carbon peaks cannot be 
detected by EDS/SEM. 



